Dynamical correlations in half-metallic manganites 
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Motivated by the recent optical and photoemission measurements on half-metallic ferromagnetic 
three-dimensional manganites, we combine a tight-binding fit of the one-particle bandstructure with 
the dynamical mean-field theory, which treats the dynamical orbital correlations and the combina- 
tion of Jahn- Teller and doping-induced disorder on the same footing. We show how all of the above 
effects are necessary to obtain good semiquantitative agreement with experimental features. As 
applications, we show how modest external magnetic fields drive drastic changes in the optical spec- 
trum, demonstrating the colossal ac magnetoconductivity. The photoemission lineshape contribution 
is evaluated, and good agreement with published experimental work is found. 



PACS numbers: 71.28-(-d,71.30-f-h,72.10-d 

In view of its potentially large technological appli- 
cations, the colossal magnetoresistance (CMR) has at- 
tracted renewed attention in the past few years. The sen- 
sitivity of transport and optical response to small mag- 
netic fields promises to open up new, attractive appli- 
cations in optoelectronic devices Q. A proper analy- 
sis of this problem requires combination of realistic ma- 
terial aspects (captured in one-electron bandstructure) 
with strong correlation effects, which are known to de- 
termine the physical response of transition metal com- 
pounds. Such approaches, recently been applied ||^ to 
the study of photoemission spectra of materials believed 
to undergo filling-driven Mott transitions, have shown 
the importance of including proper bandstructure effects 
into a state-of-the-art correlation calculation. However, 
generically, one has to deal with the added complication 
of orbital degeneracy and associated Jahn- Teller distor- 
tions, which maybe static or dynamic ||^. This is hard in 
practice, requiring one to deal with many bands-in this 
situation, a simpler approach to model the actual band- 
structure, keeping the electronically active states intact, 
is an attractive option. This is the main motivation of 
the tight-binding (TB) fits to the actual complex band- 
structure, a strategy that has been used with success for 
the case of cuprates @] . 

In this work, we use this strategy to study the 
field-dependent optical conductivity of well-doped man- 
ganites. We concentrate on the well-doped, half- 
metallic ferromagnetic state in the cubic manganite 
Lao.-rBao.sMnOs [||. We combine the TB-fit with dy- 
namical mean-field theory (DMFT) to treat the dynam- 
ical effects of local correlations with essential bandstruc- 
ture features in a self-consistent way. We show how all 
the above effects are crucial to obtain a consistent pic- 
ture of the field-dependent optical response and explicitly 
demonstrate the phenomenon of colossal optical magne- 
toconductivity. We also compute the photoemission line- 
shapes, pointing out its possible relevance to the study 
of the role of orbital correlations in the CMR materials. 

We start by noticing that in the full linear augmented 



plane wave (LAPW) Q calculation, the majority Cg 
bands, providing the conduction states, split from the 
t2g bands, but overlap with the O p bands. However, 
the two Cg bands are clearly identified along the sym- 
metry directions, making a mapping to an effective two- 
band model possible. Constructing Mn-centered Wannier 
functions with strong O pa character, only the dda = 
and ddS — t^r overlaps are nonzero. Using Slater-Koster 
tables for the cubic perovskite structure, the Cg bands 
are: = Ck ± -Dk, where ek = {t^ - ts){cx + Cy + c^), 



iter — ts) \/ C^ ~\- Cy ~\- C^ ~ CxCy 



CyC, 



and Cq = cos{kaa). The full Cg bandwidth is 2D = 
6(i(j -|- ts). To fix parameters above, notice that the dis- 
persion along r — X is equal to At^. Since the overlap 
of the majority Cg bands is omitted in the two-band fit, 
one uses the average value of the T — X dispersion, giv- 
ing Ats = -0.12, or ts = -0.03 eV. From D = 2.15 eV, 
one has = —0.69 eV. The Fermi level is set to zero 
by choosing e = 0.81 eV, completing the simple TB fit 
to the actual LAPW bandstructure. The corresponding 
total DOS is clearly different from the model DOS 
used in usual model hamiltonian treatments, and de- 
scribes the doubly-degenerate Cg bandstructure in the 
cubic perovskite geometry. Starting from similar form 
for the bandstructure Weifie et. al. concentrate on the 
polaronic physics while ignoring the undoubtedly strong 
local, orbital correlations in the FM phase It is pre- 
cisely our aim to show how a combination of strong cor- 
relations with Jahn- Teller distortion-induced orbital dis- 
order and proper one-electron bandstructure is required 
for a qualitative (in our case, we found even semiquanti- 
tative) agreement with published experimental results. 

With the single-particle dispersion relation from the 
TB-fit, the one-electron part of the hamiltonian is. 



where a = a,b label the orbital indices for the doubly de- 
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generate eg sector for the case of the cubic manganites. 
The usual spin index is dropped since one of the spin 
species is projected out of the problem in the double- 
exchange limit 1^. In this situation, the interaction part 
of the hamiltonian containing inter-orbital coulomb in- 
teraction and the Jahn- Teller coupling is given by , 

Hint = Uab ^ riialMb - ^((3^2^/ + QiS^f ) 
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The Jh ^ oo limit introduces a DE projection 
factor 1^ into the hopping term, which depends on 
the nearest neighbor core-spin correlation function: 
tij — > tij\/l-f < Si.Sj > /2S'2, and is temperature and 
magnetic- field dependent. The one-particle part is now 
written as, 

i/o= E tl\M){clc,p + h.c). 

<ij>af} 

The total hamiltonian of the system in the FM phase 
is thus H — Hq + Hint- To simplify matters fur- 
ther, we make a local rotation in the two-dimensional 
Q2 — Q3 space, so that the local quantization axis is 
parallel to z. One uses a local unitary transformation 
U^[Qi.Ti]Ui — QiTl and simultaneously transforms the 
electronic operators — UiCia- The hamiltonian 

then reads, 

^ = - E iM)UlU, {ala,,. + h.c) + Uab E "^T"a 

<ij>aa' i 
- 5 E Qi'^'^i'^ 

omitting the purely phononic part, and relabelling a^i = 
tti and ail ~ ^i- With this, we have simplified the two- 
band model in the DE limit to a Hubbard-like model 
with a disordered "magnetic field" whose source is the 
disordered Jahn- Teller distortions which exist in the well- 
doped FM state at low T. 

A plausible parameter range for the problem at hand 
is Uab > D and gQ = v < Uab, requiring one to treat 
the dynamical effects of strong local correlations and re- 
peated scattering produced by moderately strong, local 
disorder on the same footing. This is achieved by us- 
ing the DMFT as an approximation to our Sd problem. 
In d — 00, the full lattice problem is mapped onto a 
single site problem, with the "impurity" embedded self- 
consistently in a dynamical bath. In the multiband case, 
and without symmetry breaking (orbital ordering in our 
case) , the Green function and the purely local self-energy 
are functions of frequency only: Gnn'i^) = G{uj) 5 nn' and 
S„„'(a;) = S(u;)^„„'. In this situation, the DMFT self- 
consistency condition becomes, 



The remaining problem is to compute the single- 
particle self-energy, S(a;,M), in a situation where strong 
electronic correlations and disorder-induced strong scat- 
tering are simultaneously dominant, as above. We treat 
the strong orbital correlations at the local level using 
the iterated perturbation theory (IPT) away from half- 
filling, and generalized to finite temperature [lO| ] , and the 
repeated scattering effects of the local Jahn- Teller (dis- 
ordered) distortion by a proper combination of the IPT 
with the coherent-potential approximation (CPA), which 
solves the static disorder problem exactly in d — cx) [ pT| . 
In this procedure, the IPT self-energy is fed into the CPA 
local potential, which now becomes Vi = Vi — Tunti'^) — 
S(7Pa(w). The IPT local propagator is then used along- 
with this in the disorder-averaged T-matrix (CPA) equa- 
tion (but with S(ijj) — I],;„t(ijj)-|-I](7P^(a;)): the local GF 
used in the CPA equation is 



G(c.,A^) = lE 



Ei„t(w, M) - Ecpa(w, M) - el 



Solution of the CPA equation {Tii[E{Ld)]) = 0, gives 
the total self-energy corrected simultaneously for scat- 
tering caused by interactions (IPT) and disorder (CPA). 
This is fed back into the IPT subroutine that calculates a 
new self-energy, E"^™(ci;) and a new local GF, G"^™(ti'). 
These are fed back into the modified CPA routine, and 
the procedure is iterated to selfconsistency. The inter- 
acting density of states (DOS) is then obtained from the 
usual equation, p{uj,M) — —ImG{uj, M)/Tr. 

We now present the results. Below, we concentrate on 
a non-half-filled band n = 0.6 at low temperature regime 
T = O.OIL'. At "half-filling", the true state is a Neel- 
ordered orbital antiferromagnet with additional corre- 
lated distortions coming from the JT coupling. To study 
the well-doped FM regime with "melted" orbital order, 
we consider only the para-orbital state. In Fig. |l| (a), 
we show results for the local spectral function (DOS) for 
Uab ~ 2.0D. For w = 0, it shows features associated with 
the development of a correlated Fermi liquid metal as- 
sociated with collective screening of "orbital" moments, 
in analogy with what happens in the usual d = oo Hub- 
bard model. However, we observe a rich structure of the 
DOS directly related to the use of the TB-bandstructure 
instead of idealized model input bandstructure used in 
model hamiltonian treatments. 

The effect of moderate local JT disorder is modeled by 
a binary disorder distribution: P{Qi) — x5{Qi) -f (1 — 
x)6{Qi — Q), where x is the concentration of the Mn'^'^ 
sites upon doping. This is justified because hole doping 
creates locally JT-inactive Mn^~^ sites, which are ran- 
domly distributed in the host lattice of JT active M n^'^ 
ions. To make close contact with experiment we 
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choose X = 0.3. In Fig. |i| (a) we show our results for the 
two values of ty, i. e., for tij = D and Uj = \plD and for 
w = 0, J7ab/2. In particular, beyond a critical v — Vc, the 
spectrum is completely incoherent, due to strong (reso- 
nant) scattering off the disordered JT potentials. The 
DOS is characterized by a low-energy pseudogap, result- 
ing in an anomalous response. 



Bethe-Salpeter equation for the conductivity so that, 

UJ 



0.5 



S ' 
'a. 











r 


\i \ \ \ 

' \ ' A 















v=0 H,^=0 
v=UJ2 H,^=0 

v=VJ2 H^^=h^^ 









(b) 














/ /'\ \ 

u\ 




^. 
//■' 
/■•w 
/f' / 






-2 








3 



-2-10 1 

FIG. 1. (a) shows the spectral DOS for n = 0.6 and 
Uah ~ 2-D. (b) shows the corresponding photoemission line- 
shape in the non-half-filled regime for T = O.OID. The two 
values of the JT disorder and the external magnetic field are 
shown on the top of the figure. 



Next, 



we turn to the photoemission lineshape: 



/(w — fj,)p{uj). With V = Uabf^ and ti 



we see only a small "coherent" weight at the Fermi level 
(/i), as shown in Fig. |l] (b). The overall details of the cal- 
culated spectrum resemble the experimental result p^ ] 
quite well. An external magnetic field, H^xt ~ hgat re- 
sults in increased coherence, as shown by the increased 
coherent component at = ^. The role of dynamical or- 
bital correlations and JT distortions is again clear. In the 
pure DE models, one expects a large quasiparticle contri- 
bution at low-T, and the DE model combined with dis- 
ordered JT distortions yields a pseudogapped spectrum 
with no coherent contribution at low energy. Within our 
approach, we have succeeded in obtaining the quasicoher- 
ent low-energy and the incoherent high-energy satellite 
features in semiquantitative agreement with experimen- 
tal work . Notice the additional shoulder structure in 
our calculations on the higher-energy side of the coherent 
feature. This is a direct consequence of the realistic fea- 
tures in the TB-fit spectral DOS and would be missed by 
model bandstructures. Finally, in actual practice, surface 
effects and possible inhomogeneities in the sample corre- 
spond to a larger Uab, further enhancing the pseudogap 
feature in agreement with experiment. 

The optical conductivity is computed from the Kubo 
formula. In d = oo, the vertex part does not enter the 
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v=fl° Hr=A„, 

v^UJ2 H]^,^h]^„ 






v=0 




''=l7./2 


\ 

\ 




\ 

\ 

\ 

\ 

\ \ 


(b) - 


\ \ 
\ \ 
\ \ 









0.4 3 

I 



FIG. 2. (a) shows the optical conductivity for the same 
parameters as in Fig. |l| (b). (b) shows our results for the 
collossal optical magnetoconductivity for two values of the 
JT disorder. 

The calculated axx{'^) in d = oo, clearly shows up 
the magnetic-field induced spectral weight transfer in the 
"bad metal" state above *^ for x = 0.3. In Fig. | (a), 
we show Gxxi^i M) with the above parameters for two 
values of t^, which correspond to H^xt = {D{M) = D) 
and Hext = hsat {D{M) ~ \/2D), as explained above. 
With Hext = and v = Uab/2, the Drude-like part in 
axx{u}) at T > Tf^^ is strongly suppressed. This sup- 
pression is caused by the combined effects of the decrease 
of D{M) at high-T, the T-dependent reduction of the 
spectral weight (due to the T-dependence of the quasico- 
herent collective Kondo peak in our Hubbard-like model 
in d = oo) and strong scattering off disordered (local) 
JT distortions, which by itself would tend to open up 
a dip feature in the DOS. With H^xt = ^sat: we see 
a strong enhancement of the optical spectral weight in 
the low- and mid-infrared part of the spectrum. This 
is understood in terms of the physics of the d = oo 
Hubbard-like model, where Uat/ D{M) controls the dis- 
tribution of spectral weight; for large Uab/ D{M), most of 
the weight is concentrated in the high-energy part, and a 
small fraction comprises the coherent part (in the para- 
orbital metallic state), in agreement with the fact that 
the resistivity magnitude places the material in the "bad 
metal" class. With increasing D{M), high-energy spec- 
tral weight is transferred over energies on the scale of Uat 
to the quasicoherent part. This transfer scales with the 
form of D{M) as D'^{M) = D^{1 + M^)/2, completely 
consistent with experiment. 
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What is interesting is that at low energies, the frac- 
tional increase lS.axx{^, / <^xx{^,Q) is a few hundred 
percent! (see Fig. || (b)). We call this colossal optical 
magnetoconductivity (COM). Amazingly, exactly such a 
phenomenon has been observed by Boris et al. for 
Lao.eT'S'ro.aaMnOa. In practice, Hext is a few teslas, 
and, alongwith the fully spin-polarized nature of the fer- 
romagnetic metallic state, this phenomenon makes for 
interesting field-dependent optical applications. We em- 
phasize the crucial role of dynamical orbital correlations 
and JT effects here: the T and M-dependent transfer of 
spectral weight has no analogy in a non-interacting sys- 
tem. A similar calculation of axx{^) by Millis et al. 
for the case Uab = gives a completely incoherent opti- 
cal response, in contradiction with experiment, and ap- 
proaches based on the pure DE model in the cubic geom- 
etry of the perovskites jlj] overestimate the quasicoher- 
ent part vis-a-vis experiment, and completely miss the 
mid-IR features |l7| ]. 

In our picture, the extreme sensitivity of the ac re- 
sponse to a field of a few teslas arises via the change 
in the hopping, (M), which transfers optical spectral 
weight over large energy scales of 0{D). The observation 
of Simpson et al. also has a natural interpretation in 
terms of Hubbard model like physics in the orbital sec- 
tor, as shown above. In Fig. || (b), we also show how the 
COM is reduced by static JT disorder, demonstrating 
that the COM is related to the increased field-induced 
itinerance of the eg holes via tij' (M). 

In conclusion, we have shown how a combination of 
essential bandstructure aspects of the cubic perovskite 
structure with a DMFT treatment of strong orbital cor- 
relations in the eg sector in the DE limit describes the 
optical and photoemission spectra of the metallic CMR 
manganites in a semiquantitative way. Additional inter- 
esting applications of our modelling to investigate field- 
dependent magneto-optical response across the para- 
ferro transition are being studied, and will be reported 
elsewhere. Our treatment can be extended to include 
other broken symmetries in the spin and orbital sectors, 
and should be generally applicable, with suitable modifi- 
cations, to other half-metallic TM-oxide ferromagnets. 
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